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Abstract 

Bogoliubov's theory states that self-interaction effects in Bose-Einstein condensates produce a 
characteristic linear dispersion at low momenta. One of the curious features of Bogoliubov's theory 
is that the new quasiparticles in the system are linear combinations of creation and destruction 
operators of the bosons. In exciton-polariton condensates, this gives the possibility of directly ob- 
serving the negative branch of the Bogoliubov dispersion in the photoluminescence (PL) emission. 
Here we theoretically examine the PL spectra of exciton-polariton condensates taking into account 
of reservoir effects. At sufficiently high excitation densities, the negative dispersion becomes visi- 
ble. We also discuss the possibility for relaxation oscillations to occur under conditions of strong 
reservoir coupling. This is found to give a secondary mechanism for making the negative branch 
visible. 
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Recent experimental advances have achieved the condensation of exciton-polaritons in 
semiconductor microcavity structures [l-3]. The short hfetime of the exciton-polaritons 
on the order of picoseconds means that the condensate is rather different in nature to a 
traditional atomic Bose-Einstein condensate. The condensate exists only resulting from the 
replenishing of the polaritons from a reservoir of uncondensed polaritons, which is in turn 
populated by illumination by a laser. Despite this difference, such condensates exhibit many 
of the characteristics expected in equilibrium condensates, ranging from superfluidity ji] to 
vortex formation . In the work of Ref. l7|, it was shown that self- interaction effects of 
the condensate cause the dispersion characteristics of exciton-polariton condensates follow a 



Bogoliubov dispersion relation. Although in the work of Ref. j?] no negative branch of the 
Bogoliubov dispersion js] was detected, recently a four- wave mixing experiment has revealed 
the presence of the negative branch 9j. The four- wave mixing experiment was originally 
theoretically proposed in Ref. 10||. The relative difficulty of the observation of the negative 
branch was attributed in this work to the bright condensate emission which easily masks 
the weaker negative branch. 

In this paper we present a detailed theoretical analysis of the photoluminescence (PL) of 
the negative Bogoliubov branch. In contrast to the four- wave mixing approach of Ref. jlO|, 
the PL is calculated directly via two-time correlation functions of the polariton equations of 
motion. In particular, we incorporate the effect of the bottleneck polaritons which is known 
to strongly influence the dispersion of the polaritons HI]. To this end, we first reformulate 



the theory of Ref. 



Ill in a Heisenberg-Langevin formalism. The reformulation makes it 



clear that the theory is a modified Bogoliubov theory defining new bosonic excitations in 
the system. 

The polaritons are assumed to obey the Hamiltonian 
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where p'^{x) is a creation operator for a polariton at position a?, nji{x) is the number operator 
of the reservoir polaritons, m is the mass of a polariton, g is the interaction energy of the 
polaritons, and g is the polariton-reservoir interaction. For example, including only the 
exchange interaction g = ^^^^ |X|^^, where e is the electronic charge, is the Bohr 
radius, X is the excitonic Hopfield coefficient, and A is the sample area, and we use SI units 



12| . The polaritons then obey the Heisenberg-Langevin equations of motion 



^h^^ = [pix), H]+'- iRiuRix)) - ^)p{x) (2) 



where R{nR{x)) is the stimulated gain coefficient [13| of the reservoir polaritons into the 
jjh is the decay rate of the polaritons through the microcavity mirrors. 



condensate anc 

Following Ref. 11|, we linearize equation ([2]) such that it only involves terms linear in either 
the polariton or reservoir operators. To achieve this goal, we expand the reservoir operator 
into its average and fluctuation components 

nnix) = 4 + ^6n{x) (3) 

Wo 

where n^j is the average reservoir number and ipo is the amplitude of polaritons in the 
condensate. The prefactor of 6n{x) is chosen for later convenience. Substituting ([3]) into 
([1]) and rewriting the operators in terms of their Fourier components p{x) = 'YlikPk^^^'^ 
and 5n{x) = ^ J^k^^k^^'"^ ^ obtain 



k,k' ,q 



k,q 

We note here that 6nk should be interpreted as the amount of density fluctuations in the 
reservoir with Fourier component k, and should not be confused with the number of reser- 
voir polaritons with momentum k. We now follow the same procedure to the Bogoliubov 
prescription to obtain a Hamiltonian that is bilinear in the variables (pfc,p^fc, Suk) by picking 
out terms in the summation which involve polariton operators with k = 0, and set these to 
their average values Pk=o — ^ = V'oc"*^^*'^'* jl^, where /ly = ^n5j + 5'|\E'oP is the condensate 



energy. This gives 



H =||v&or + E [(^ + 9< + 29\^o\')pipk 
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The resulting equation of motion for the polariton operators is thus 

= + 9nl + 2g\mA p^ + g^lpl^ 
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dt \ 2m 
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Meanwhile, the reservoir equation of motion obeys 

'^''^^^^ - G{x) - ^nriRix) - R{nn{x))p^{x)p{x), (7) 



dt 



where G is the Langevin noise operator for the number operator for reservoir polaritons [15 1 
and 7_r//i is the decay rate of the reservoir polaritons. The noise operator G originates from 
the coupling of the reservoir modes to high energy excitations induced by the laser pump. 
We have neglected the diffusion of the reservoir polaritons since they have a negligible effect 
on the dispersion characteristics. Substituting into (^^, expanding in Fourier space, and 
performing the Bogoliubov linearization we obtain 



^Suk = ^{Gq- ^Rn\)5{k = 0) - 7ij(5nfc 
dt njj 



- ^^^^^(^opI + ^SP-.) - i?'(4)|^op5r.„ (8) 

where we have assumed a homogeneous pump G{x) = Gq. 

The k ^ components of eqns. (|6]) and ([8]) can be conveniently summarized in the form 

3 

dt 



where = {pk^P^k^ ^^^^^^Z^, (5nfce ^t^rt/n^^ ^^^^ 
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The matrix Mk is identical to that given in Ref. 

In Fig. [1] we show the real part of the eigenvalues of Mk for typical parameter values. 
We identify two regimes which depend primarily on the relative magnitude of the scattering 



rate R and the reservoir decay rate 7ij. In Ref. [ll|], it was generally assumed that '^r was 
large, giving the characteristic fiat dispersion regime in the vicinity of = 0, as seen in 
ig. [1^. However, considering that reservoir polaritons originate from the bottleneck region 



16l | which generally have a longer lifetime than the condensate polaritons, it is possible to 
consider the opposite regime where 7^ is small and the scattering rate R is large. In such a 
regime (see Fig. [T)d) the eigenspectrum shows a split dispersion at A; = 0, with a rounding of 
the Bogoliubov spectrum, such that the particles once again acquire an effective mass. This 




FIG. 1: (Color online) (a) (b) Energy eigenvalues and (c) (d) ratio of coefficients I'^^fc P/l^^^ P for 
the eigensolutions of the matrix (solid lines). Dashed lines show the corresponding values for 
Bogoliubov theory (no reservoir), with energy dispersion Esogik) / Eq = \J {kaY{{kaY + ~^^) and 
coefficients — 1/2 = + 1/2 = (^^ + g^l) /2EBog{k). The energy scale is measured 

in units of Eq = i^^, where a is the experimental length scale (e.g. in GaAs a = lO^'^cm and 
£"0 = 0.68meV). Parameters used are gmP^ = 1, gijjQ = 1, R'{inP^inP^ = 1, R{n^) = 1, g = g, and 
(a) (c) = 1 (corresponding to a flat dispersion regime) (b) (d) 7^^ = 0.1 (corresponding to a 
relaxation oscillation regime). All parameters in units of Eq. 

is caused by relaxation oscillations in the system, where the reservoir and the condensate 
repeatedly exchange population if displaced out of equilibrium. A simplified description can 
be obtained by expanding the Gross-Pitaevskii and reservoir equations in Ref. [Ill] around 
their steady state values ip{t) = ipo + Si/j{t) and riRit) = + {n'j^/ilJo)6n(t) and ignoring 
interaction effects g = g = 0. For small k we obtain 

The imaginary part of the eigenvalues gives the oscillation frequency, which is 

u = \^R'{n'n)R{nlWo ' ilR + R'{n'M?/^- (10) 

In the limit of small scattering R{nii) or large 7^^ the frequency becomes pure imaginary, 
indicating that only damping occurs. 
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The eigenvalues of the matrix give the new collective excitations of the system. 
Specifically, new quasiparticle operators 

feM = E%v. (11) 

may be defined where are the eigenvectors of My.. The matrix allows for two types 
of solutions, either dispersionless or bosonic, defined as whether the commutation relation 

K.^u] = \<?-W^? (12) 

is zero or non-zero respectively. Here we used the identities = 5n_y. and [dn^., = 0. 

Operators with non-zero commutators may be normalized to ±1, defining new bosonic modes 
of the system. From Fig. ^ we see that no bosonic solutions exist in the fiat dispersion 
regime, since all eigenvectors satisfy Im^P = l^^fc^P, which corresponds to [fefc,^^^] = from 
( !T2|) . In the the regime beyond the fiat dispersion there is always one dispersionless solu- 
tion with P = l^ife^P corresponding to a renormalized density fluctuation solution. The 
solution with Im^P > Im^^P corresponds to a solution with [&fc,&fe^] = 1, which is a new 
"dissipative" Bogoliubov destruction operator. The solution with jw^^p < \u'^^\^ meanwhile 
corresponds to [6^^,6_fe] = —1, the dissipative Bogoliubov creation operator (see also Fig. 
[T]i). Putting in the correct time dependences, outside the fiat dispersion regime we may 
associate the solutions to be 6^ = (6^, fol^e"^*'"^*/'', ^mfcC"*'"^*/''), corresponding to the pos- 
itive, negative, and dispersionless branches respectively. Since ( ITTl) appears to admix both 
bosonic and number operators one may wonder how to interpret such an operator. The 
clearest interpretation is that this is a displaced Bogoliubov operator = b'f^ + u]^6nk, 
where the first two components form the boson operator and the reservoir creates displace- 
ments in the vacuum from this state. In this case the amount of reservoir fluctuations of 
wavelength k displaces the vacuum for Bogoliubov particles b'f^ at momentum k. 



To calculate the PL spectrum we apply the methods presented in Ref. 12|]. The PL 
spectrum is given by 

PL{k,u) oc \Ck\^Re[{pi{t = 0)pfc(a;))] 

= n\'Re[Y,{uir<' {bl\t = 0)K{oo))] (13) 

where Uk is the inverse of Uk (normalized to satisfy bosonic commutation relations), Ck is 
the Hopfield coefficient for the photonic component of the polaritons, and tildes denote time 



Fourier transformed variables. In the diagonal basis the equations of motion (Q are 

'^iF = EiJ^ + (14) 

where we have added a Langevin noise operator = (F^, 0) to account for a thermal 
population of dissipative Bogoliubov particles |14i], and are the eigenvalues of M^. The 



noise operator is assumed to obey correlations of the form 13 1 



{blFl) = 5aTknt (15) 

where n^^ = l/(e^'=/'=sT _ ^ Re[El] = -Re[El], = lra[El] = lm[El], T is the 

temperature, and 6ij is the Kronecker delta. This ensures that the dissipative Bogoliubov 
particles obey thermal statistics {h'^jfik) = n)^, with other off-diagonal correlations vanishing. 
The choice of ( |T5i) is chosen primarily because it is the simplest choice, and is possible to 



generalize to other distributions 17|. 



After a Fourier transform of (HM and inserting the correlations we obtain the expression 
for the PL 



PL{k,u) oc ICfcpRe 



(16) 



hjjj — — iTk hu + — iT^ 

Here the first term corresponds to the positive dispersion, weighted by the thermal pop- 
ulation and the dissipative Bogoliubov coefficient, while the second term is the negative 
dispersion, which only appears when there is appreciable Bogoliubov mixing between the pk 
and p'_i^ operators. In Fig. [2^ and[2]D we plot the PL spectrum corresponding to low and high 
density regimes respectively. At low density we see that only the positive branch is visible. 
Here there is negligible mixing between the components (ITT]) due to the small off-diagonal 
components of M^. The positive branch is populated via the thermal noise field Fk, while the 
negative branch remains dark. At sufficiently high density, the off-diagonal components of 
Mk become large enough such that there is some mixing of the components ( ITTj) . This results 
in a visible negative dispersion in the PL spectrum. A simple criterion for when the negative 
branch becomes visible may be derived. As may be seen by inspection of the Bogoliubov 
expression for the factor Im^^P = + g^o) f^Esogik) — 1/2, the negative component is 
only appreciable for momenta < gijQ. Due to the finite linewidth of the dispersion, the 
negative dispersion only becomes resolvable beyond momenta > Tk- Thus the negative 
dispersion only becomes visible for densities exceeding the criterion gi/j^ ~ F^. 



7 




12 3 4 
ka 



1 


2 3 


4 




ka 






1 2 3 4 1 2 3 4 

ka ka 



FIG. 2: Photoluminescence of the excitation spectrum of exciton-polariton condensates for various 
parameters. Regimes of (a) just above condensation threshold (b) thermahy depleted high density 
(c) high density and low temperature (d) zero interactions are shown. Parameters used are gn^ = 1, 
R'{n\)n\ = 1, /?(nOj) = 1, 5 = ^, 7^ = 1, and (a) = Qmin% keT = 2, (b) V'g = n%, keT = 2, 
(c) V'o = '^fi' ^bT = 0.1, (d) g = g = 0, TpQ = n^, ksT = 2. A chemical potential of /x = — 1 was 
assumed in the thermal distribution in order to account for finite size effects of the condensate. All 
parameters in units of Eq. Zero detuning and a Rabi splitting of lOE'o was assumed to calculate 
the Hopfield coefficient Ck- 

It is interesting that for low temperatures where n)^ is small, only the negative excitation 
branch is visible in the PL spectrum (Fig. [2b). The reason for this can be understood by 
the following argument. According to the inverse relation of flTTl) . the loss of a polariton out 
of the system (which is the basic process underlying the PL) is described by the operator 

Pk = ui'bk + u^^hl^e-'"'^'^'"' + v'ibmue-'^^'l^ (17) 

The energy change of the system associated with the first term is /iT + Cfe, which is nothing 
but the standard mechanism for the PL with a positive dispersion. With no thermal pop- 
ulation of dissipative Bogoliubov particles, the first term automatically gives zero, thus the 
positive dispersion remains dark. The second term is associated with the gain of a dissipative 
Bogoliubov particle and the loss of two condensate particles, which has an energy change of 
— (/iT + Cfc) = /iT — Cfc- Unlike the first term where a dissipative Bogoliubov particle 
needs to be originally present, the second term does not require this condition and occurs 
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ess of the initial population. We note that a similar effect has been observed in atom 



18| . The third term causes a loss of a condensate particle with an energy change of 



fiT independent of k, giving a dispersionless spectrum. This can in principle give rise to 
a fiat PL emission as seen in Fig. [T] if there are strong enough density fluctuations in the 
reservoir, although we do not assume that this is typically the case in current experimental 
systems. 

There is another mechanism for making the negative branch visible, which is in a regime 
where reservoir scattering is strong. To see this effect, we show in Fig. [2jl the PL spectrum 
for the illustrative case of zero interactions g = g = 0. We again see that the negative 
branch is visible along the eigenvalues of M^. In this case the negative branch is visible 
not because of self-interactions, but due to mediation via the reservoir mode. From our 
simulation results we flnd that the negative branch becomes visible for the same condition 
as that given to observe relaxation oscillation (i.e. that (fTOj) is real). The mechanism for 
this is that due to the relaxation oscillations there is an effective coupling of the pk and 
plfc operators mediated by the 6nk operator. In practice both the self-interactions and the 
reservoir effect are likely to contribute to making the negative dispersion visible. 

In summary, we have calculated the PL spectrum of the excitations of exciton-polariton 
condensates. At sufficiently high densities such that gipQ ~ Ffc, the negative branch of the 
Bogoliubov spectrum becomes visible. The negative dispersion should be visible even if no 
thermal population is present, due to the nature of the PL emission being associated with 
the loss of a polariton from the system. In the regime of long reservoir decay times, a regime 
of relaxation oscillations was identified, which was also found to be able to illuminate the 
negative dispersoin. One question which remains is why a four-wave mixing experiment 
was required in order to see the negative dispersion in Ref. |9l], instead of spontaneously. 
One explanation is that the bright emission from the condensate itself may be masking the 
negative branch, and only the faint tail at high momenta can be accessed 10|]. Another pos- 
sibility is the reduction of the condensate self-interactions at higher density due to reduction 



of the Bohr radius 



19 



20| . This may contribute to make the negative branch weaker than 



expected due to the suppression of Bogoliubov mixing. 
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